Frozen Quasi-Long-Range Order in the Random Anisotropy Heisenberg Magnet 
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Extensive Monte Carlo simulations are used to investigate the low-temperature properties of the 
random anisotropy Heisenberg model, which describes the magnetic behavior of amorphous rare- 
earth-transition metal alloy. We show that the low-temperature phase in weak-anisotropy region is 
characterized by a frozen-in power-law spin-spin correlation. Numerical observation of the power-law 
exponent rj indicates non-universal behavior. 
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Note: 

Owing to the page limitation, some parts of this manuscript 
are omitted in the published version. To obtain the published 
version, see the comments at the beginning of the TeX file. 



Power-law correlation, or quasi-long-range order 
(QLRO), usually emerges just at the critical tempera- 
ture of phase transition. The notable exception is the 
Kosterlitz-Thouless phase of two-dimensional XY model, 
in which power-law correlation persists down to zero 
temperature^ Recently it has been found that the low- 
temperature QLRO phase is quite common in disordered 
three-dimensional systems^ such as the so-called "Bragg- 
glass phase" of impure superconductors^*! nematic phase 
of liquid crystal in a porous mediaj£*£ and amorphous 
rare-earth-transition metal alloyi These systems are de- 
scribed either by random-field or random-anisotropy spin 
models. Harris, Plischke and Zuckermann first used the 
following random anisotropy model to study amorphous 
metal magnelA 
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where Si is a Heisenberg spin on the lattice site i of 
simple cubic lattice, J > is a ferromagnetic coupling 
constant, D > is the strength of uniaxial anisotropy, 
the former sum runs over all nearest-neighbor pairs, and 
rii is a random unit vector which describes the random 
anisotropy at the site i. Similar model, in which Si ■ Sj is 
replaced by (Si ■ Sj) 2 , is used to study disordered liquid 
crystals. While an Imry-Ma-type argument of model 
iJTJl leads to a conclusion that any nonzero anisotropy 
D destroys long-range magnetic order in three dimen- 
sions, it is predicted by field theoretical analysi s 2 ! 10 ^ 1 
that the model realizes QLRO at low-temperature and 
weak-anisotropy region. The phase transition between 
this QLRO phase and paramagnetic phase is studied by e- 
expansion techniques,- 2 and it has been found that there 
are no renormalization group fixed point in a replica- 
transformed model of usually this is interpreted as 
a sign of discontinuous transition, but no experiments 
report it. Numerical studies of model have been 



restricted to approximate discretized models or small 
and the QLRO phase was numerically con- 
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firmed in Rcf. 15 . Similar phase was also found in the 
simulation of disordered liquid crystal^ The possibility 
of spin-glass transition in the strong-anisotropy region 
has been investigated by Monte Carlo simulations, but 
no conclusive result has been obtainedii&i^ Experimental 
studies have shown that there is a "spin freezing transi- 
tion" below which zero-field cooled and field cooled mag- 
netization differiiSiii This phase is referred to as "cor- 
related spin glass" phase, 18 owing to the developed spin 
correlation length as opposed to the usual spin glasses. 
Similar glasslikc behavior was found in the recent exper- 
iment of disordered liquid crystal^ However, the pre- 
dicted QLRO has not yet been reported experimentally. 20 

In the present work, we numerically investigate 
the low-temperature properties of model on an 
L x L x L lattice for several values of D/J (including 
infinity) using the exchange Monte Carlo methodiSi 
Details of the simulation are summarized in Table 
|U We simultaneously simulate N r identical systems, 
assigning different temperature to each one. At each 
Monte Carlo step, exchange of states between adjacent 
temperatures Tj and T,+i is tried and accepted with 
a probability min[f ,exp {{E t - E i+l ){l/T t - 1/T i+1 ))], 
where Ei denotes the energy of the state at Tj. 
Thus each state random-walks the temperature 
space. To maximize the diffusion constant of this 
process, a set of temperature {T{\ is adjusted so that 
{{E i+1 )-(Ei))(l/T i+1 -l/Ti) ~ -I. The average energy 
can be easily estimated beforehand from simulations 
of smaller systems. Equilibrium is checked by dividing 
the measured data into groups and calculating averages 
for each group, then discarding the data of initial 
groups which are not equilibrated. In the Metropolis 
update of large D/J cases, the new spin is chosen with 
a probability proportional to exp(D/T(ni ■ Si) 2 ) and 
acceptance probability is calculated using only the first 
term in ^ For small D/J cases, over-relaxation type 
update is also employed When measuring the spin- 
glass order parameter, three replicas are placed at each 
temperature and overlaps between them are observed. 
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Note: 

In the previous versions of this manuscript, the following 
observation scheme was used in some of exchange MC 
simulation. This scheme is incorrect, since the marked 
replica had always been at higher temperature just be- 
fore it is observed, and average over such samples is bi- 
ased to high-temperature side. 

A special observation scheme is used when we investi- 
gate a large system at temperature far below the critical 
point: a replica is "marked" when it visits the highest 
temperature, and if a marked replica reaches the lowest 
temperature, we store its configuration into memory and 
"unmark" the replica. After MC steps proportional to 
iVjj, a large number of independent, equilibrium configu- 
rations at low temperature are obtained. 



Average over random anisotropy configurations is per- 
formed over 48 < N s < 900 configurations. Statistical 
errors are estimated from sample-to-sample fluctuations 
by the Jackknife procedure. All the simulations were per- 
formed on Fujitsu VPP5000 vector processors at JAERI. 
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TABLE I: Details of the exchange Monte Carlo simulations. 
Nr, MCSj, MCSo, and Ns denote number of replicas placed 
between T m i n and T max , number of Monte Carlo steps used 
for thermalization and discarded, number of Monte Carlo 
steps used for observations, and number of samples used for 
averaging over randomness, respectively. The sign "*" refers 
to all the values of D/J used in Fig. 5. 

We observe both magnetization m = L~ 3 J2i &i an d 
overlap parameter q = L~ 3 1h where q\= S" ■ Sf is 
an overlap at site i between two independent replicas a 
and (3 which we run in parallel. Figure ^ shows plots of 
Binder's cumulant of fh and q: 
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for the D / J = 4.0 case, where (■ • ■) and [• • ■] denote ther- 
mal and sample averages, respectively. A clear crossing 
at T w 1.423 can be observed for both B m and B q which 




FIG. 1: Plots of B m and B q for the D/J = 4.0 case. 



clearly indicates that there is some kind of order at low 
temperatures (Note that the transition temperature is 
only 2% smaller than that of D — case, as was found in 
Ref. — ). To investigate the nature of the ordered phase, 
we observe the effective exponent of rj defined between 
two different system sizes L\ and L2: 



2-d 



M(m(L 2 ) 2 )]-ln[(m(£ 1 ) 2 )] 
7hS ~ ln(L 2 )-ln(ii) ' ( ) 



where d denotes spatial dimension (in the present case, 
d = 3). When spin-spin correlation at a distance r is pro- 
portional to r 2 ~ d ~ r> , then r] e g asymptotically approaches 
to 77 as the size is increased, while when the correlation 
is short ranged, [(to 2 )] is proportional to L~ d , therefore 
i] c ff — 2 is approached. Figure [21 shows plots of rj g for 
the D/J~ 4.0 case: The value of rj c g is almost indepen- 
dent of temperature below the critical point. This be- 
havior is quite different from that of Kosterlitz-Thouless 
(KT) phase^i in which rj e ff jumps from 2 to 1/4 at the 
critical point, then gradually decreases down to as the 
temperature is lowered down to T = 0. The difference 
from the KT phase becomes remarkably clear when one 
sees a probability distribution of the overlap parameter 
P(q') = [(S(q -q'))]. The inset of Fig. shows P(q) 
of two sizes L = 8 and 12, for the D/J = 4.0 case, at 
a very low temperature T / J = 0.6815 which is well be- 
low the critical point; the sharp peaks, which become 
sharper as the size is increased, indicate that only one 
kind of state is dominant in the ordered phase and there 
are no critical thermal fluctuations. From these results, 
one can depict the nature of the low-temperature phase 
as follows: the system is trapped around one of the two 
ground states ±{S ? | ''}, in which spin-spin correlation de- 
cays as Si -Sj ~ \i— j'| _7?_1 . At the transition tempera- 
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FIG. 2: Plots of the effective exponent r/ e ff against tempera- 
ture for the D/J — 4.0 case. Statistical errors are smaller than 
the size of the symbols. The inset shows probability distribu- 
tion of the overlap parameter q at D/J = 4.0, T/J = 0.6815. 



ture, the system "freezes" into this QLRO state, then the 
magnetic QLRO and the spin-glass order simultaneously 
emerge. Baker and Kadanoff suggested similar QLRO 
phase for models with finite zero-point entropy, such as 
antiferromagnetic Potts models, in which thermal fluctu- 
ation is present even at T = 0. They speculated that in 
these models the point T = 0, along with the whole low- 
temperature phase, is renormalized to some nontrivial 
fixed point. But in the present case, the behavior of P(q) 
suggests that there is only one kind of low-temperature 
state and thermal fluctuation is irrelevant in the low- 
temperature phase. 

Next, we consider the low-temperature phase in the 
strong- anisotropy region, in which QLRO is expected to 
be destroyediiii^ Figure |21 and 21 shows plots of B m 
and B q for D/J = 10 and D/J = oo cases, respec- 
tively.The case D/J = oo correspond to a kind of Ising 
spin-glass model Y^tij) Jij a i a j with Jy = Hi ■ ftj. The 
plots of B m shown in Fig. |3| and 0] do not show 
crossing behavior and are consistent with the past re- 
sults , and indicate that the spin-spin correlation is al- 
ready short ranged at D/J = 10. As for the spin- 
glass order, however, the crossing of the plots of B q is 
rather ambiguous; they seem to collapse into a line in 
the low-temperature region. This result is suggestive of 
Kosterlitz-Thouless-like quasi-long-ranged spin-glass or- 
der, but it should be remembered that numerical simula- 
tions of the three-dimensional ± J spin glasses could not 
exclude a (false) possibility of KT-like transition, unless 
recently large-scale simulations have become possible. 23 
So, in the present work we reserve conclusion on this is- 
sue. The phase diagram shown in Fig. [S] summarizes 
these results. In the weak-anisotropy limit, magnetic or- 
der is either slowly decaying QLRO or long ranged, which 
is difficult to distinguish numerically. 

Now let us consider the following question: Is the ex- 
ponent 7] universal all over the QLRO phase? The field 
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FIG. 3: Plots of B m and B q for the D/J = 10.0 case. 




FIG. 4: Plots of B m and B q for the D/J — oo case. 



theoretical 4— e expansion analysis 211 suggests that a sta- 
ble fixed point at D/J = 0(e) governs the nature of the 
whole QLRO phase. Figure Eta) schematically depicts 
the finite-size behavior of ry e // expected from this theo- 
retical framework. Although the 4 — e expansion analy- 
sis only foretells the presence of a stable RG fixed point, 
there should be another, an unstable fixed point since the 
QLRO is destroyed in the strong-anisotropy limit. When 
the anisotropy ratio D/J is smaller than a threshold (un- 
stable fixed point), it is attracted to a stable fixed point 
as it is renormalized, while it will diverge when larger 
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FIG. 5: Phase diagram of the random anisotropy magnet, 
suggested by the results of the present work. Dotted lines 
and shaded area are ambiguous phase boundaries which are 
hard to identify by numerical simulations. 





FIG. 6: (a) Schematic representation of finite-size behavior 
of the effective exponent r? e ff expected from renormalization- 
group analysis. Thin and thick curves are finite-size values 
of smaller and larger systems, respectively. The bold straight 
lines show the infinite-volume limit, (b) Effective exponent 
rjefi of several values of D/J, measured at T/J = 0.5 for all 
cases. The data represented by a dotted line are calculated 
from correlation function of spin configurations obtained by 
simulated annealing of L — 64 system. 



than the threshold. Thus plots of rj e ff against D/J 
at a low temperature for different sizes will intersect at 
two fixed points. Figure EJL) shows plots of ?7 e // mea- 
sured at a very low temperature T/J = 0.5 (far below the 
transition line T ~ 1.4) for several values of D/J and for 
different sizes. In a region 3.0 < D/J < 6.0, the two 
plots coincide within statistical errors, unlike Fig. EJa)- 
This suggests a nonuniversal QLRO phase in which the 
asymptotic value of r\ continuously changes with varying 
D/J, or at least that the renormalization-group flow is 
very slow. 

To see if the observed nonuniversal QLRO is asymp- 
totic or not, we investigate larger system L = 64 with 
simulated annealing (SA) method^ The apparent draw- 
back of SA is that it does not guarantee equilibration, 
therefore obtained configuration may not be the true 
ground state. Energy of the obtained states is given in 
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TABLE II: Energy per spin of the annealed states, averaged 
over randomness. The third row shows standard deviation of 
sample-to-sample fluctuations. 



Using the configurations obtained by SA, we cal- 
culate spin correlation function defined as G(r) = 
(3L 3 ) -1 J2\i-j\= r &i ' Sji where the summation runs over 
all pairs whose relative position is either (r, 0, 0), (0, r, 0), 
or (0, 0, r) lattice spacings. For each random anisotropy 
configuration, only one annealed configuration is gener- 
ated and G(r) is averaged over ten random anisotropy 
configurations. The summation over lattice sites signif- 
icantly reduces statistical errors, and this small number 
of random anisotropy configurations is enough to obtain 
reasonable precision. Figure shows plots of G(r) for 
various values of D/J. For the cases of D/J < 4.0, 
G(r) in a range 4 < r < 16 is well fitted by a power- 
law form G(r) = ar - * 7-1 , while for D/J > 5.0 a form 



G(r) 



= n.r-V-1 



exp(— x/£) is well fitted. Owing to the 



periodic boundary condition, the data in a region r > 16 
deviate from the fitting forms. The values of r\ for each 
D/J obtained by the fitting are shown in Fig. ??, and 
they coincide with the previous results within statistical 
errors. The inset of Fig. [7]shows plot of rj and l/£ against 
D/J: G(r) becomes short ranged at around D/J = 5.0, 
therefore the QLRO phase should persist at least up to 
this anisotropy strength. 

Now that we have found a QLRO phase at around 
D/J — 4.0, which is much larger than that of ex- 
perimentally realized materials D/J ~ 1.0, a question 
arises: Why has this QLRO phase never been observed 
in the past experiments, such as small angle neutron 
scattering2£? One possible reason may be the temper- 
ature dependence of the anisotropy D: Generally D in- 
creases as the temperature is lowered, while J remains 
almost constant. 18 Since the ground state is completely 
modified when D/J changes, reaching the true ground 
state requires complete reorganization of spins, which is 
unlikely to happen no matter how slowly the system is 
cooled. 

In summary, we have carried out extensive Monte 
Carlo simulations of the three-dimensional random 
anisotropy Heisenberg model and found a QLRO phase 
which is characterized by frozen power-law spin corre- 
lations in the low-temperature, weak-anisotropy region. 
This phase persists at least up to the anisotropy strength 
D/J w 5.0. Finite-size behavior of the effective power- 
law exponent in this phase indicates nonuniversal behav- 
ior, contrary to the field theoretical prediction of univer- 
sal exponent. 
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FIG. 7: Spin-spin correlation G(r) in the annealed states for 
D/J = 3.0, 4.0, 5.0, 6.0, 8.0, and 10.0 cases (from top to 
bottom), averaged over ten samples for each case. Lines are 
results of fittings, either G(r) — ar - * 7-1 (solid lines) or G(r) = 
ar~ n ~ 1 exp(— r/f) (dotted lines). The inset shows plot of the 
fitting parameter r\ and l/£ against D/J. 
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